Collective excitations in trapped boson-fermion mixtures: from demixing to collapse 
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We calculate the spectrum of low-lying collective excitations in a gaseous cloud formed by a 
Bose-Einstein condensate and a spin-polarized Fermi gas over a range of the boson-fermion cou- 
pling strength extending from strongly repulsive to strongly attractive. Increasing boson-fermion 
repulsions drive the system towards spatial separation of its components ("demixing"), whereas 
boson-fermion attractions drive it towards implosion ("collapse"). The dynamics of the system is 
treated in the experimentally relevant collisionless regime by means of a Random-Phase approxima- 
tion and the behavior of a mesoscopic cloud under isotropic harmonic confinement is contrasted with 
that of a macroscopic mixture at given average particle densities. In the latter case the locations 
of both the demixing and the collapse phase transitions are sharply defined by the same stability 
condition, which is determined by the softening of an eigenmode of either fermionic or bosonic 
origin. In contrast, the transitions to either demixing or collapse in a mesoscopic cloud at fixed 
confinement and particle numbers are spread out over a range of boson-fermion coupling strength, 
^"v and some initial decrease of the frequencies of a set of collective modes is followed by hardening as 

evidenced by blue shifts of most eigenmodes. The spectral hardening can serve as a signal of the 
impending transition and is most evident when the number of bosons in the cloud is relatively large. 
We propose physical interpretations for these dynamical behaviors with the help of suitably defined 
CO ' partial compressibilities for the gaseous cloud under confinement. 

CN 

PACS numbers: 03.75.-b, 67.60.-g 



CO 

o 
o 

CN 



c3 



^ , I. INTRODUCTION 

in 

Dilute boson-fermion mixtures are currently being produced and studied in several experiments by trapping and 
cooling vapors of mixed alkali-atom isotopes [II. 0- [J 01I3- El The large variety of combinations of atomic species and of 
choices of magnetic sublevels, together with the additional possibility of tuning the interactions by means of Feshbach 
resonances provide a unique opportunity of investigating extensively the effects of boson-fermion interactions, both 
in the case of mutual repulsions and in the case of attractions. 

The boson-fermion coupling strongly affects the equilibrium properties of the mixture and can lead to quantum 
phase transitions. Boson-fermion repulsions can induce spatial demixing of the two components when the interaction 
*^ ■ energy overcomes the kinetic and confinement energies 0, and several configurations with different topology are 
possible for a demixed cloud inside a harmonic trap |S|. Although spatial demixing has not yet been experimentally 
observed, the 6 Li- 7 Li experiments of Schreck et al. [l| appear to be quite close to the onset of the demixed state 
. In the case of attractive boson-fermion coupling the mixture becomes unstable against collapse when the boson- 
, induced fermion-fermion attractions overcome the Pauli pressure 0, El ■ Collapse has been experimentally observed 
in a 40 K- 87 Rb mixture 0. Fermion pairing into a superfluid state has also been predicted in the case of attractive 
interactions [TTj . but for spin-polarized fermions this is expected to happen in the p-wave channel at temperatures 
much lower than the current experimental limit. This possibility will not be considered in the present work. 

The dynamical properties of boson-fermion mixtures have been previously investigated theoretically mainly in the 
mixed phase, both for the macroscopic homogeneous system and in mesoscopic clouds inside a harmonic 

trap. In the latter case these studies have used a sum-rule approach 0,E1, perturbation theory ^t|> or a Random- 
Phase approximation (RPA) 0, ^ • 

In this paper we investigate the effect of boson-fermion interactions on the spectrum of collective excitations as 
the mixture ap pro aches an instability and look for the dynamical signatures of the approaching instabilities. In a 
previous paper |21| we have studied the transition to demixing of a confined cloud in the collisional regime and found 
an analytical condition for the dynamical transition point. We focus here on the experimentally more relevant case 
of the collisionless regime, where we adopt an RPA scheme. Within this formalism it is possible to study at the same 
level the evolution of the cloud towards demixing on one side and collapse on the other, thus providing a unified 
picture of very different states. We characterize the dynamical properties of the mixture not only by its spectrum but 
also through its generalized partial compressibilities. 

We comparatively examine a homogeneous mixture and a mixture in external harmonic confinement. Whereas in 
the former case the RPA spectra are symmetric under changing the sign of the boson-fermion coupling, inside the 
trap we find different behaviors on going from the case of boson-fermion repulsions to that of attractions. This is due 
to the fact that in the theory the transition in the macroscopic homogeneous system is approached at fixed density of 
the two species, while in the trapped mesoscopic cloud we work in the realistic situation of constant particle numbers 
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and fixed external confinement. 

The paper is organized as follows. In Sec. [n] we review the conditions for demixing and collapse as obtained from 
static considerations, and in Sec. IHII we introduce the RPA formalism for the collective excitation spectrum. Sec. II VI 
describes the predictions of the RPA for the case of a homogeneous system, while Sec. reports the results for the 
collective excitation spectrum of a mixture under external harmonic confinement. Finally, Sec. lVII contains a summary 
and our main conclusions. 



II. EQUILIBRIUM STATE: CONDITIONS FOR DEMIXING AND COLLAPSE 



We consider a dilute boson-fermion mixture at zero temperature. The Hamiltonian which describes the system is 

h 2 v 2 

(1) 



PC"' V 



+ V e l t (r)- t , c 



where and *S?f are the usual bosonic and fermionic field operators, V^'f (r) are the external confining potentials, 
(J>B,F are the chemical potentials, and m B}F the atomic masses. We have introduced the couplings g BB — 4:TTh 2 a BB /m B 
for boson-boson interactions and g BF = 2nh 2 a BF /fn r for boson-fermion interactions, in terms of the s-wave scattering 
lengths cibb and clbf an d of the reduced boson-fermion mass m r . We will always assume o,bb > 0. We are here 
neglecting the fermion-fermion interactions since we are considering a spin-polarized Fermi gas where s-wave collisions 
are forbidden by the Pauli principle and higher-order collisional couplings are negligible at the temperatures of present 
interest. In the dilute limit we also neglect the quantum depletion of the condensate. 

The equilibrium state of the mixture is described at mean field level by the self-consistent solution of the Gross- 
Pitaevskii equation for the condensate extended to include the boson-fermion interactions, 



h 2 v 2 

2m B 



Kxt( T ) + 9bb n B (r) + g BF n F {r) 



and by the Hartree-Fock equations for the fermions coupled to the boson density, 

h 2 v 2 



2m F 



+ Vext( r ) +9BF n B (v) 



xjji(r) = £itpi{r) 



In Eqs. (01 and J2J the condensate density is given by 

n B (r) = |$s(r)| 2 



and the fermion density by 



n F (r) 



\^(r)\ 2 6(ti F ~e t ). 



(2) 



(3) 



(4) 



(5) 



The chemical potentials are obtained by imposing the normalization conditions Nb,f — J d 3 f it-b,f( v ) giving the 
numbers of particles of the two species. Excited-state condensate wave functions and Hartree-Fock orbitals will also 
be needed in the construction of the RPA dynamical susceptibilities in Sec. II I II 

Under isotropic harmonic confinement V^f \r) — rriB ,f^b f 7 " 2 /^i f° r the values of Np ~ 10 4 as in current exper- 
iments we have verified that the fermionic equilibrium profiles are well described by the solution of the generalized 
Thomas-Fermi approximation (TFA) 



h 2 V 2 7n7(r) ^ ,2/3 



6m F y/n F (r) 



+ An z /°{r) + V L F xt (r) + g B Fn B {r) = hp, 



(6) 



which includes the surface kinetic energy effects in the form of von Weiszacker Here A = {Ti 2 /2m F )(67r 2 ) 2 / 3 . The 
fermionic density profiles for various values of the boson-fermion coupling gsF < are shown in Fig.^ The result of 
the TFA is practically indistinguishable from the full solution of Eqs. J3J and JSJ, since the shell effects play only a 
marginal role in the spherically symmetric case for the above choice of fermion number. The inset in Fig. ^ reports 
the phase diagram of the harmonically confined mixture for the case g BF < (see below). 
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At increasing values of the boson-fermion coupling the mixture becomes unstable against demixing (in the case 
9bf > 0) or collapse (in the case g BF < 0) 0. In a homogeneous system at fixed particle densities a linear stability 
analysis based on a mean-field energy functional predicts the same condition for collapse and for demixing p^ : 

9BB9FF - 9bf = ° ( 7 ) 

— 1/3 

with g FF = (2/3) An F playing the role of an effective fermion-fermion coupling due to the Pauli pressure in the 
Fermi gas. 

In the presence of external harmonic confinement the effects of finiteness substantially modify the description of 
the instabilities. In the case of demixing the transition is smoothed out and the onset of demixing with increasing 
boson-fermion coupling is signalled first by a diminishing boson-fermion overlap energy 0. The condition of such 
partial demixing is given by 

9bb \ N 2 J 5 NT 



where the parameters are a = (15 3 ^ 5 /A8 1 / 2 )(m F /m B ) 3 ^ 2 (a BB /a ho ) 3 ^ and c 2 = (48 1 / 3 /15 3 / 5 )(6/7r) 2 / 3 (a BS /a /lo ) 2 / 5 , 
with ciho = (Ji/msLUs) 1 ^ 2 ■ On further increase of the boson-fermion coupling one reaches the point where the fermionic 
density vanishes at the center of the trap, which occurs at 

dyn 

£b£L = ^L. (9) 

9BB fJ-B 

This condition also corresponds to a sharp upturn in the frequencies of collective excitations in the collisional regime 
|2l| . If the boson-fermion coupling is further increased the overlap between the two clouds becomes negligible, the point 
of full spatial separation being well predicted by Eq. J7J) in a local-density approximation . At full separation the 
most likely configuration with current experimental parameters e.g. in a 6 Li- 7 Li mixture Q is the "egg" configuration, 
formed by a core of bosons surrounded by a shell of fermions. 

In the case of attractive interactions, at increasing values of the boson-fermion coupling the cloud is predicted to 
enter first a region of metastability where the zero-point oscillations still prevent collapse. For equal masses and 
trapping frequencies the condition of metastability as derived by Miyakawa et al. ^(J is given by 

g BB -2^N F [W > 

On further increase of \gBF | the mixture reaches the point of collapse, that in the numerical solution of the equilibrium 
equations J2J and © we localize at the point where convergence fails. We have verified that this instability point is 
well described by Eq. Q in a local-density approximation Q, E| j 

\9bf 1 \_\ ™ b 1 1 1/2 

9BB 



_2m F (67r 2 n jF (0)) 1 /3 asB _ 

where n F (0) is the fermion density at the center of the trap. From the conditions illOfl and it can be seen that 
for large boson numbers N B ^S> Np the metastability region is absent and the cloud is predicted to go directly from 
a stable to an unstable configuration. The inset in Fig. ^ shows the phase diagram of the confined mixture in the 
case g BF < with the choice of parameters that we shall adopt later in the study of the RPA spectrum of collective 
excitations. Similar curves for the critical boson number needed to reach collapse have been reported by Roth p3l |. 

In summary, the criterion obtained in Eq. Q for a macroscopic cloud is still useful to locate within a local-density 
approximation the end of the transition to either demixing or collapse in a boson-fermion mixture inside a harmonic 
trap. However, the transitions in a mesoscopic cloud kept under fixed external confinement are spread out and their 
approach will best be characterized by dynamical signatures, as we proceed to show below by spectral calculations in 
cases of increasing boson-fermion coupling. 



III. THE RANDOM-PHASE APPROXIMATION FOR A BOSON-FERMION MIXTURE 



The RPA is equivalent to a time-dependent Hartree-Fock theory of the dynamics of a quantum fluid in the col- 
lisionless regime |24| and was first introduced for dilute boson-fermion mixtures under confinement by Capuzzi and 
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Hernandez |T^. It treats on the same footing the density fluctuations of the two atomic species, satisfying the /-sum 
rules and allowing for Landau damping of the bosonic excitations by the Fermi cloud. Since it neglects correlations 
between density fluctuations, the RPA is valid in the dilute limit specified by the inequalities nsa%g <C 1 and 
&f|o_bf| ^ 1; where kp = (6n 2 nF ) 1 ^ 3 is the Fermi wave number. 

The RPA yields the spectrum of collective excitations in the linear regime from a set of coupled equations for the 
density fluctuations Suf and Sub, which are obtained by assuming that the fluid responds as an ideal gas to external 
perturbing fields 8Up and 5Ub plus the Hartree-Fock fluctuations. The RPA equations in Fourier transform with 
respect to the time variable read 

Sn F (v,Lj)= J d 3 r' X QF (r,r',Lu)[6U F {r',L;)+g BF 5n B {r',Lu)} (12) 

and 

5n B (v,uj) = [ d 3 r'x° B (r,r',u)l5U B (r',Lu)+g BB Sn B (r',L;)+g BF Sn F (r',Lu)]. (13) 



In the inhomogeneous system consistency between equilibrium and dynamics requires to take into account static 
mean- field interactions. Thus x° F m Eq. I|12l) is the Lindhard density-density response function constructed with the 
Hartree-Fock orbitals from Eq. ©, 

X 0F (r,r» = £ /( g i) ~ /fe) ^( r )^(r0^ W|(r) (14) 

with /(e) = 6(hf — £) and r] = + . Similarly, x° B m Eq. (|13|) is the condensate density-density response function 
constructed with the excited states <pi of the Gross-Pitaevskii equation with energy eigenvalues Ei , 



X 0B (r,r» = £ 



$ B (r)g| (r')0* (r% (r) d> B (rQg| (r)# (r)& (rQ 
lu - (Ei - n B ) + ir) lu + (Ei - fx B ) + ir/ 



(15) 



The RPA Eqs. i|12|l and (|13|) are most easily solved in terms of a matrix of density-density response functions, which 
are defined as 

Xaa . (r, r', t - = JP4t4 = - t')([5p*(r, t),5 Pcr ,(r', t')}), (16) 

where Sp a (r, t) are the density fluctuations operators in the Heisenberg representation. Upon introducing the matrices 



XFF XFB \ -o _ f X° F 



and G = ( ^ ^ BF ) the RPA equations read 



XBF XBB J ' X V X° B J ' " V 9BF 9BB 

*(r,r»=x°(r,r» + J d 3 r" X °(r, r", w) G X (r", r', lu). (17) 
In fact, the boson-boson interactions can be resummed into the Bogoliubov response functions x s ° 9 , defined from 

X Bog (r, r', U ) = x° B (r, r', u) + J <?r" X 0S (r, r", w ) x s °V, r', w). (18) 
Thus the RPA equation can alternatively be expressed as 

X(r,r»=x°(r,r»+ / d 3 r" X °(r, r", to) G X (r", r', w) (19) 



with x° = ( X n Boo ) and G = ( ^ '^? F I . In this way the excitations of the bosonic cloud are described as 



x B ° 9 J "~\9bf 
a gas of Bogoliubov phonons interacting only with the fermions. 
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IV. SPECTRA OF THE HOMOGENEOUS MIXTURE 



The RPA equations for the homogeneous fluid mixture with average particle densities ns and rip can be solved 
algebraically in Fourier transform with respect to the relative position coordinate r — r'. The solution takes the form 

1 / X° F (Ku) g B FX Bog (k,uj)x 0F (Kuj) 



where the "dielectric function" e(k, uj) is the common denominator of the four partial responses and is given by 

e(k, ui) = 1 - g% F x B ° 9 (K uj) x° F (K w). (21) 

The spectrum of collective excitations is obtained by searching for the poles of the response matrix or equivalently by 
imposing e(k, ui) — in the complex frequency plane. In Eqs. I|20|l and l|21ll the Bogoliubov response function is 

U) (ui + irj) — c B k z — [nk z /2niB) 



with cb — {^f)BB n B I vhb)^^ i and the Lindhard response function is 

X 



0F (k ui] = V ^ ep) - f( 6p+ ^ (23) 
4^ w - (e P +k - e p ) + irj 



where £k = hk 2 /2mp. The summation in Eq. I|23|) can be carried out explicitly (see e.g. the book of Fetter and 
Walecka |25|) to yield the spectrum of single particle- hole pair excitations and the corresponding reversible part of 
the susceptibility of the ideal Fermi gas. 

Equations l|2U|) and i|21|) were first introduced by Yip [lj , who evaluated numerically the spectra in the case of weak 
boson-fermion coupling. For cb > up, where vf = hkp/mF is the Fermi velocity, the spectrum at gsF = contains 
a stable bosonic phonon and a fermionic particle-hole continuum. When the boson-fermion interactions are turned 
on these modes hybridize and repel each other, as can be analytically shown in the long wavelength limit. From the 
condition e(k, u>) = 0, using Eq. (|23) and x° F (k,u;) -► n F k 2 /(m F oj 2 ) for lo > VFk and k — * 0, we obtain the correction 
to the Bogoliubov sound velocity due to boson-fermion interactions as c = cb[1 + g B F m BnF/(g%B m F' n B)] 1 ^ 2 ■ This 
result agrees with that obtained by second-order perturbation theory by Viverit and Giorgini |l4| . In the case cb < vf 
the bosonic sound falls instead inside the fermionic particle-hole continuum and its broadening by Landau damping 
must therefore be considered. 

At increasing values of the boson-fermion coupling one must resort to a numerical solution of the RPA equations. 
Figure |21 shows the fermionic spectrum Imx_FF(k, ui) for various values of gBF, the taking other parameters from the 
values at the center of the trap in a 6 Li- 7 Li mixture with equal trap frequencies u>p = ^>b = ui a = 2tt x 1000 s -1 
and bosonic scattering length asB = 0.27 nm. While in the RPA all four partial response functions Xaa' (k, uj) share 
the same poles, we have chosen this spectrum because it shows most clearly the role of the fermionic particle-hole 
continuum. 

We consider first the case cb > vp- As the gas approaches the instability points for demixing or collapse as given 
by Eq. ||7J|, we observe a softening of the particle-hole spectrum until most of the oscillator strength is transferred 
to the bosonic mode at the point of collapse or demixing (see the first panel in Figure^. In the case cb i$ vp we 
observe that the boson-fermion interactions help to stabilize the bosonic mode, since mode repulsion pushes it out of 
the particle- hole continuum (see the second panel in Figure This is not possible however in the case cb <§C Vp, 
where the effect of boson-fermion hybridization is to break the particle-hole continuum into two parts (see the third 
panel in Figure yj. In this case as the system approaches the instability point it is the original bosonic mode which 
softens and broadens. 

The softening of the spectrum is directly reflected in the partial "compressibilities" , obtained from the long wave- 
length limit of the static susceptibilities as 

tv n m r du' lm X ^(k^0,uj') 
Xr7 a>{k^0,uj = 0)= . (24) 



In the homogeneous mixture these quantities can be calculated analytically from the limiting behaviors of x Bog and 



X 0F : 



1 



X(k^ 0,a; = 0) = = —I ~ 9BB 9BF ). (25) 

gBBgpF - g B F V SBF gFF 
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From the above equation we immediately see that the RPA predicts a divergence in the compressibilities at the same 
instability points as obtained from static considerations. Alternatively, the instability point of the homogeneous 
mixture can be predicted from the emergence of unstable modes according to the condition 

lim e(k, iw) = (26) 

at long wavelengths. 

Finally, from the form of the RPA solution in Eqs. I|2()(l and l|21|l it is natural to define the effective fermion-fermion 
and boson-boson coupling as gf/f = g 2 BF x B ° s (k, uj) and gfM = g\ F x° F (k, u>), respectively |2^,|23|. From this point 
of view the instabilities occur when the effective interactions, which are negative in the limit (k — ► 0, u> = 0), overcome 
the Fermi pressure for the fermions and the boson-boson repulsions for the bosons. 



V. SPECTRA OF A MESOSCOPIC CLOUD IN A TRAP 



We consider now a mixture under external harmonic confinement, choosing for simplicity the spherically symmetric 
case V^f = niB.FOJQr 2 /2 with the same trap frequency u>o for the two atomic species. As for the homogeneous gas, 
we investigate how the spectrum of collective excitations is modified as the confined cloud moves toward demixing or 
collapse. Important differences arise from the presence of the confinement, under which both the Bogoliubov sound 
and the fermionic particle-hole spectrum are quantized. In particular the discrete energy levels depend on the static 
Hartree-Fock fields and shift as we vary the coupling strength, so that we know a priori that the spectra will not be 
symmetric under a change in the sign of gsF- 

For the finite system we follow the usual convention of evaluating the spectral strength functions Xo-o-' ( w ) [l8l,l28| . 
which are defined as 

WH = fd 3 r ( d 3 r' 5U;,(r)x™'(r,r> ' ,u)5U a (v') (27) 



where SU a (r) are the intensities of the time-dependent external potentials driving the two atomic species. In the 
case of spherical confinement these have the form 8U&(r) = fa{r)Yi m (f) and correspondingly we can decompose the 
response functions into components of definite angular momentum, 

4-7T 

W(r,r» = £ — x ^,( r y,^ m (f)Y4(f'). (28) 

l,m 

In these equations Yi m (f) are the spherical harmonic functions and r denotes the unit vector along r. 

Phase separation into the "egg" configuration or collapse has symmetry corresponding to I = and therefore we 
focus on the monopolar excitations of the type SU a (r) — Uq (r/aho) 2 ■ This type of perturbation can easily be produced 
in current experiments since it corresponds to just a modulation of the trap frequencies. In this case the strength 
function takes the form 

x!T» = ^ / d 3 r [ d 3 r'^-^ X l ZHry,u;). (29) 

J J a ho a ho 

We have obtained these RPA spectra by numerically solving the integro-differential problem posed by Eqs. I|12(l - I|15|l 
in conjunction with the determination of excited-state orbitals for the condensate and for the Fermi gas. A few points 
on the numerical solution are worth noting: (?) the orbitals ipi and fa are discretized on a uniform mesh of about 
250 points expanding to approximately twice the size of the density profiles; (ii) convergence is achieved both in 
the chemical potentials and in the respective energy spectrum; and (Hi) for attractive boson-fermion coupling the 
numerical solution is computationally more demanding since a larger basis set is needed to properly describe the 
equilibrium density profiles as we approach the collapse. 



A. Across demixing 

We have calculated the RPA spectra of a mixture of 7 Li and 6 Li atoms Q with trap frequency luq = 27r x 1000 s -1 , 
bosonic scattering length agg = 0.27 nm, and a tunable clbf chosen to span the phase diagram from the mixed 
state to partially demixed states (as we shall see below, the spectrum shows a dynamical transition well before full 
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separation). Experimentally the value of the boson-fermion scattering length can be tuned e.g. by using a Feshbach 
resonance. 

In Fig. El we show the imaginary part lmx l pp(u) of the strength function for Np = 10 4 and various values of Nb- 
In the absence of boson-fermion interactions the cloud responds to a monopolar excitation with only two well-defined 
modes, a purely fermionic one at frequency u>f — 2 u>q |29| and a purely bosonic one at frequency uj\, ~ \fh u>o |3C|. The 
latter, of course, is not present in the fermionic spectrum for cibf = 0. Upon turning on the boson-fermion interactions 
the bosonic peak appears in the fermionic strength function with a considerable oscillator strength, showing that due 
to the boson-fermion coupling the fermionic cloud can oscillate in resonance with the bosons even in the collisionless 
regime. This "bosonic" mode moves to higher frequency as clbf is increased, analogously to what is found in the 
homogeneous mixture. A second "bosonic" peak also appears in the fermionic strength function, corresponding to a 
mode at lo^i — \/i4wo in the absence of boson-fermion coupling. This n = 2 mode is not excited by a monopolar 
perturbation in a pure bosonic cloud, but is excited in the mixture due to the dynamical coupling between bosons 
and fcrmions. Within the RPA the fermionic density fluctuations act as further perturbing fields on the bosons, and 
vice versa. . 

For values of the frequency close to the original purely fermionic peak we observe a strong fragmentation of the 
strength function, most of which is already present at mean-field level in the Hartree-Fock fermionic spectrum. 
This fragmentation, which has also been reported for other choices of system parameters [T^ . [l9| , is due to boson- 
fermion interactions lifting the angular degeneracy of the single-particle levels for fcrmions in the trap. In addition 
to fragmentation we observe a softening of the original fermionic spectrum. The softening stops when the topology 
of the cloud changes, with the opening of a hole in the fermionic density at the trap center, although strong overlap 
of the two clouds is still present. This dynamical transition has also been found in the collisional regime (2l| and 
occurs in both regimes at the same point, which is well described by the analytical expression @. For larger values 
of the boson-fermion coupling the original fermionic component of the spectrum moves upward, again in analogy 
with the behavior of a trapped mixture in the collisional regime [2 it . This effect is intrinsically due to the harmonic 
confinement: at increasing asF the fermion cloud forms a shell around the bosons with increasing average density 
and decreasing thickness, thus increasing its stiffness in response to the monopolar perturbation. In the homogeneous 
system instead the softening of the fermionic spectrum continues up to the point of full thermodynamic separation. 
The role of the number of bosons in the cloud is also evident from Figure |21 

As in Sec. If VI it is possible to monitor the softening or hardening of a part of the spectrum in the trapped cloud by 
looking at the "minus-one" moment of the monopolar spectral strength distribution. This spectral moment is specially 
sensitive to the low-frequency part of the spectrum and, owing to the Kramers-Kronig relations, yields the static limit 
X 1 ^®' ( w — 0) of the strength function, thus defining four generalized compressibilities which measure the reaction of the 
cloud to a radial pressure exerted on a given component. In Fig. 0] (left panels) we show the fermionic compressibility 
X l pp{uJ = 0) both for positive and negative values of gsF and for three choices of the number of bosons. For positive 
qbf we observe an initial increase of Xpp{u> ~ 0), followed by a slow decrease after the dynamical transition has 
taken place. The position of the maximum corresponds to the dynamical point of demixing in Eq. (J5J). 

At variance from the homogeneous case, the four compressibilities in the trapped cloud have a nonsymmetric behav- 
ior for positive and negative values of the boson-fermion coupling. In particular this means that the compressibilities 
cannot be predicted by a simple local-density approach. For example, in the homogeneous gas one always finds a 
negative value of ~xbf{^ = 0) for gsF > irrespectively of the actual values of the densities, whereas as is shown in 
Fig. 0| (right panels) this compressibility becomes positive in the confined cloud after the dynamical demixing point. 

B. Towards collapse 

We consider now the case of attractive boson-fermion interactions, choosing the trap parameters of the 6 Li- 7 Li 
mixture as in the previous section. As is shown in the inset in Fig. ^ when the numbers of bosons and fermions are 
comparable the cloud can be found in a metastable state, while for larger numbers of bosons a wide stable region is 
present in the phase diagram. The collapse is reached at smaller values of cibf as the number of bosons is increased, 
as can be understood from the fact that the boson-mediated fermion- fermion attraction increases with Nb- 

In the homogeneous gas the RPA spectra for attractive boson-fermion interactions coincide with those given in Fig. [21 
for the case of repulsive interactions. This does not apply to the trapped mixture, as can already be seen from the 
compressibilities in Fig. Contrary to the behavior of the homogeneous gas where the generalized compressibilities 
diverge at the instability points (see Eq. I|25| l). in the harmonic trap for negative values of clbf they decrease and for 
large numbers of bosons tend to vanish as the cloud approaches collapse. This is an indication that spectral weight 
is being transferred to high-frequency modes, and the effect becomes more dramatic as the number of bosons in the 
mixture is increased. In the limit Nb 3> Np the decrease of the fermion-fermion compressibility can be analytically 
predicted using a single- mode approximation (see the inset in Fig.^J. The details are given in Appendix lAl 
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A blue shift of the modes in the fermionic spectrum is found already at the Hartree-Fock level, as is illustrated in 
Fig. |S]by reporting the density p%{s) of excited pairs, 



at increasing Nb for Nf = 10 4 and various values of cibf- In Eq. 1)30(1 we have chosen only the pairs with A^ = to 
obtain a monopolar particle-hole excitation spectrum, with An = 1 since this gives the most important contribution 
at low frequencies. 

The spectra resulting from a full RPA calculation of the monopolar fermion-fcrmion response in the region of 
boson- fermion attractions are shown in Figs. and {7\ When the numbers of bosons and fermions are comparable we 
observe a shift of the pure fermionic peak towards lower frequencies. We interpret this shift as being associated with 
the lifting of the degeneracy of the single-particle Hartree-Fock spectrum by boson-fermion attractions, which allow 
transitions between states that are closer in energy. At increasing values of Nb we find instead that the fermionic 
mode is blue-shifted by a much larger amount, as a result of the static Hartree field of the bosons. The amount of the 
shift can be predicted analytically in the limit Nb S> Nf, where the interactions with the bosons simply renormalize 
the trap frequency (see Appendix^). In the same limit the bosonic modes are almost unaffected by the interactions 
with the fermions, as the fermionic Hartree field becomes almost negligible for the bosons. 



In summary, in this paper we have studied how the mutual interactions of bosons and fermions affect the spectrum 
of collective excitations of a dilute boson-fermion mixture in the collisionless regime as the mixture goes towards 
either demixing or collapse. To this purpose we have used the Random-Phase approximation, which is most suited in 
the dilute regime since it accounts correctly for the limit of small couplings and generally conserves the /-sum rule. 
We have studied both a homogeneous mixture and a mixture in isotropic harmonic confinement, where new effects 
arise due to the static mean-field interactions. We have indeed found that both the boson-fermion coupling and the 
confinement strongly change the spectra of collective excitations. Our main results are as follows. 

First of all, in the presence of boson-fermion interactions the excitation of one component of the mixture will set into 
motion also the other component. There also is hybridization of the excitation modes, e.g. the fermionic component 
can also oscillate in resonance with the bosonic cloud at the frequency of a bosonic mode and vice versa. 

In the homogeneous gas we observe qualitatively different behaviors depending on the ratio between the Bogoliubov 
sound velocity cb and the Fermi velocity vp- In the case cb > Vp at increasing values of the boson-fermion coupling 
we find mode repulsion between the broad particle-hole continuum and the Bogoliubov sound, and the fermionic 
spectrum softens as the gas approaches either instability point. When cb is just below vf the Bogoliubov sound is 
damped at low couplings, but a "revival" of this mode can take place as the boson-fermion interactions are increased. 
Softening and damping of the bosonic mode is instead observed when cb <C vf- Within the RPA we have also derived 
a dynamical condition for the two instabilities, which coincides with that obtained by studying the linear stability of 
the mean- field energy functional. 

In the trapped system we have studied the I = monopolar excitations. While in the absence of interactions the 
spectrum of density fluctuations shows only two sharp and independent modes, the boson-fermion coupling induces 
fragmentation of the original fermionic mode, revealing the presence of many occupied single-particle levels in the 
Fermi sphere. At increased repulsive couplings we observe a softening of the spectrum up to the point where the 
topology of the fermionic cloud changes by the formation of a hole at its center. The spectrum is then shifted to 
higher values of the frequencies at larger values of cibf- For negative values of the coupling we observe some frequency 
decrease for comparable numbers of bosons and fermions, but as the number of bosons increases the static Hartree 
field of the bosons drives a sizable blue shift of the fermionic modes. 

From the comparative analysis of the collective excitation spectra of the mixture in the homogeneous state and 
inside a trap we have learnt that the inhomogeneous trap potential can change quite drastically the physical picture, 
since the instabilities are approached along thermodynamically different routes. In the homogeneous gas one works 
at fixed average particle densities, while the confined cloud must be treated at fixed external confinement. We have 
found that the partial compressibilities cannot be described by a local-density picture, and while they are expected 
to diverge in the homogeneous gas on approaching the instability points, in the trap they saturate to finite values 
as the component of the mixture are spatially separating and vanish as collapse is approached. This suggests that 
a description of boson-induced fermion-fermion interactions, e.g. in the study of the instability against pairing, may 
also be affected by the confinement. 

In this work we have neglected the static and dynamical effects of correlations and the quantum depletion of the 
condensate These should become important in the proximity of collapse as the gas densities increase. The 
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effects of temperature may also be of interest, with the thermal depletion of the condensate and the blurring of the 
Fermi surface yielding, for instance, damping of the Bogoliubov sound. Finally, new dynamical effects may arise when 
the components of the trapped mixture have very different atomic masses and will be studied elsewhere. 
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APPENDIX A: SINGLE-POLE MODEL FOR THE TRAPPED MIXTURE 

We derive in this Appendix a single-pole model to describe the partial fermion-fermion compressibility for attractive 
boson-fermion interactions in the limit of large boson numbers. This model accounts rather well for our numerical 
results and shows how the inhomogeneous harmonic confinement in combination with the boson-fermion coupling can 
give rise to a completely different behavior in the compressibility relative to the homogeneous case. 

In the limit N B ^> N F the dynamical coupling between bosons and fermions is quite small, but the fermionic density 
is largely affected by the static mean field of the bosons. The bosonic equilibrium profile is only weakly affected by the 
presence of the fermions, so that in the Thomas Fermi approximation we can write ns(r) = [/xs— Vj^. t (r)]6(B.B— r)J g B B 
with R B = (2/jb/jtibWq) 1 / 2 . In these conditions it is useful to introduce an effective potential acting on the fermions 
as 32] 

V e F ff( r ) = \ m F^lr 2 +9BFn B {r) 



m B g BF 

1 6(R B - r) 

rriFgBB 



—— ^bQ(Rb - r). (Al) 

9BB 



For the choice N B = 2.4 x 10 7 and N F — 10 4 the bosonic and fermionic clouds have the same radii, so that 
the effective potential is seen by the fermions at all r and acts as a harmonic trap with renormalized frequency 
Co = w [l - (m B g BF )/(m F g BB )} 1/2 . 

An analytical estimate of the compressibility can be obtained by making a single-pole approximation on the fermion- 
fermion strength function 

Imx FF (w) = --A5(u) - co e ) (A2) 

where for monopolar excitations we take io e — 2d). The oscillator strength A is determined from the /-sum rule, which 
for a trapped system reads |3^| 

- / — hLoImXFFiu) = Mi = / d 3 r\VSU F (r)\ 2 n F (r) (A3) 

J tt m F J 

In our case 5U F = U${r / ato) 2 and the RHS of Eq. I|A3|) is determined by the mean square radius of the equilibrium 



fermionic density: M\ = AhcuoUfi (r 2 ) / at (m B /m F ). This is readily evaluated from the Thomas-Fermi approximation 
rip(r) = (6Tr 2 )~ 1 [(2m F /h 2 )(fi F — m F Co 2 r 2 /2)} 3 / 2 9(fj, F — m F Cb 2 r 2 /2)), where for R F ~ R B the chemical potential is 
simply given by /i F = TiCb^Np) 1 ^ . We therefore obtain for the oscillator strength the value 



A = 34/32-2/3^4/3 ^V 2 ^) . (A4) 

Finally, from the imaginary part of the strength function we can calculate the compressibility by means of the 
Kramers-Kronig relation, with the result 

Xff(. = 0) = ^NTU 2 (^(^)\ (A5 ) 
27T^ v u ) \m F J 
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This depends on obj? through Co and gives a quantitative agreement with the numerical calculations (see the inset in 

Fig mi. 
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FIG. 1: Equilibrium density profiles df(r) (in units of aj^f) as functions of the radial coordinate r (in units of dho) for the 
fermionic component of a 6 Li- 7 Li mixture under isotropic harmonic confinement, as obtained from the solution of the Hartree- 
Fock equations (dots) and from the Thomas- Fermi approximation (solid line, hidden by the dots). The curves from top to 
bottom correspond to the choices cibf = —4 nm, (ibf = —2 nm, and cibf = 0. The other parameters are Nb = 10 5 , Nf = 10 4 , 
cibb = 0.27 nm and u>o = 2-7T x 1000 sec -1 . The inset reports the phase diagram of the mixture with boson-fermion attractions 
in the adimensional variables o,bf /cibb (linear scale) and Nb/Nf (log scale), showing the regions of stability (S), metastability 
(M) and collapse (C). The crosses show the points where the spectra of the trapped mixture are explored in this paper. 
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FIG. 2: Imaginary part of the fermion-fermion density response function — ImxFF(k, uj) for k = 0.01 &f as a function of u) 
(in units of Tik 2 F /iriF) in a homogeneous 6 Li- 7 Li mixture for various values for the coupling asF- the solid, dotted, dashed, 
and dash-dotted lines correspond to four equally spaced values of |osf| from zero to the demixing or collapse point lying at 
asF = 7.9 nm. The arrow marks the location of the Bogoliubov sound mode in the absence of boson-fermion interactions. The 
three panels from top to bottom correspond to the choices cb/vf ^ 1.4, cb/vf ^ 0.92 and cb/vf ^ 0.14. The corresponding 
values of the particle densities (in units of o^) are n F = 8.9 and n_g = 32555, 14649, and 351.5 from top to bottom. 
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FIG. 3: Spectrum of monopolar response (in log scale and arbitrary units) as a function of frequency to (in units of loo) for a 
boson-fermion mixture in harmonic confinement with Nf = 10 4 and Nb = 10 5 (top panel), 10 6 (middle panel), and 2.4 x 10 7 
(bottom panel), for various values of clbf (in nm) as indicated in each panel. The vertical dotted lines indicate the frequency 
of fermionic (left) and bosonic (right) modes in the absence of boson-fermion coupling. 
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FIG. 4: Generalized compressibilities (in units of Xff at obf = 0) of a boson- fermion mixture in harmonic confinement as 
functions of the boson- fermion scattering length bbf (in nm) for Nf = 10 4 and Nb = 10 5 (top), 10 6 (middle), and 2.4 x 10 7 
(bottom). The small and large arrows indicate the points of onset of partial demixing and of the dynamical transition from 
Eqs. JHJ and respectively. The lines are a guide to the eye. The inset shows an enlargement of the region with o,bf < in 
comparison with the analytical predictions from Eq. I|A50 (dashed line). 
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FIG. 5: Hartree-Fock fermionic spectral functions (in arbitrary units) as functions of lo/loq for Np = 10 4 and Nb = 10 5 (top), 
10 6 (middle), and 2.4 x 10 7 , at various values of obf as indicated in the panels. 
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FIG. 6: Spectrum of monopolar response (in log scale and arbitrary units) as a function of frequency u) (in units of u>o) for a 
boson- fermion mixture in harmonic confinement with Nf = 10 4 and Nb = 10 5 (left panel) and 10 6 (right panel), at various 
values of clbf (in nm) as indicated in each panel. The vertical dotted lines indicate the frequency of fermionic (at integer 
multiples of 2ojq) and bosonic modes in the absence of boson-fermion interactions. 




FIG. 7: Spectrum of monopolar response (in log scale and arbitrary units) as a function of frequency lu (in units of u>o) for 
a boson-fermion mixture in harmonic confinement with bbf = —2 nm and N F = 10 4 at various values of Nb- Nb = 10 
(top), Nb = 10 6 (middle) and Nb = 2.4 x 10 7 (bottom). The vertical lines indicate the frequencies of fermionic modes (dotted 
lines) and bosonic modes (dot-dashed lines) in the absence of boson-fermion interactions. For the above sets of parameters 
thermodynamic collapse is predicted to occur at ubf — —4.9 nm for Nb = 10 5 , —4.4 nm for Nb = 10 6 , and —4.0 nm for 
iV s = 2.4 x 10 7 . 



